REACTION ENHANCED DIFFUSION IN SPHERICAL MEMBRANES 
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Abstract 

The reversible reactions like A + B <=> C in the many-component diffusive system affect the diffusive properties 
of the constituents. The effective conjugation of irreversible processes of different dimensionality takes place due to the 
stationarity in the system and can lead to essential increase of the resulting diffusive fluxes. The exact equations for the 
spatial concentration profiles of the components are difficult to treat analytically. We solve approximately the equations 
for the concentration profiles of the reaction-diffusion components in the spherical geometry in the application to the 
problem of the enhanced oxygen transfer through a biological membrane and to the mathematically similar problem of 
surface diffusion in a solid body. In the latter case the spherical geometry can be an adequate tool for describing the 
surface of a real solid body which can be modeled as a fractal object formed of sequences of spherical surfaces with 
different radii. 

Keywords: Enhanced diffusion; Transport phenomena; Thermodynamics; Membranes; Reversible chemical reactions 
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1 Introduction 

The phenomenon of "facilitated through chemi- 
cal reaction diffusion" is familiar to various domains 
of science. Its essence consists in the diffusion en- 
hancement through the intermittance of some chem- 
ical reaction involving the diffusive components. A 
straightfoirward example is the so called "vacancy 
enhanced diffusion" of doped impurity into a solid 
body. Since long ago it was widely known that the 
presence of crystal imperfections is able to facilitate 
the impurity penetrating thus considerably enlarging 
the effective diffusion coefficient. From the thermo- 
dynamical point of view the facilitation comes from 
the presence of several fluxes of the diffusant (as if it 
were propagating through several different channels, 
or rather consist of several species having different 
diffusivity) and intervenience of the reversible chem- 
ical reaction between them. Thus, in a solid body 
the impurity atioms are known to reside either in the 
nodes of the crystalline lattice, (substitutional uimpu- 
rity) or in the space between nodes (interstitial impu- 
rity), the latter being fast diffusants and the former - 
slow ones. The reversible interchange of these two 
species tooks place and is facilitated by the presence 
of vacansies and eigen interstitial atoms of the ma- 
trix. An another example is taken from the biology 
domain and it concerns the phenomenon of the fa- 



cilitated oxygen transfer through cellular membranes 
Ffl . The principal reaction scheme consists in re- 
versible "tying" the ligande molecules by some slow 
macromolecular carrier. The facilitated transport of 
oxygen is possible via some fermentative kinetics 
(which should not however affect the chemical prop- 
erties of O2) and especially through reversible reac- 
tion with haemoglobine or myoglobine. In what fol- 
lows we refer for concreteness to this example though 
our consideration can be applied to the vast variety of 
situations falling into the same reaction scheme. 
The reaction process with Hb and oxygen can be (al- 
though very schematically) represented by the fol- 
lowing expression 



Q 2 + Hb^± Hb0 2 



(1) 



which means the formation of an (unstable) complex 
HbO^, the rates of forward and reverse reactions are 
k\ and 

The system of balance equations for three con- 
stituents can be written as follows: 



DAc = p + q 
D p Ac p = p 



(2) 



D p Ac c 



-P 



where c, c„ and c c stand for concentrations of O2, 



Hb and complex HbC>2 respectively; D and D p are 
diffusion coefficients for O2 and Hb (assuming that 
both Hb and #&0-complex possess the same dif- 
fusion coefficients since the Hb molecule is much 
larger in size than oxygen one and both pure Hb and 
the complex HbC>2 must have similar diffusive prop- 
erties). Instead Hb some other protein, e.g. Mb can 
appear. The consuming term q (here set to 0) should 
be added, for example, if treating the problem of mus- 
cular oxygen transport with myoglobine. The reac- 
tion rate p from elementary kinetics reads as follows 

p = k\cc p — k_c c (3) 
imposing minimal nonlinearity on the system. 



ary conditions for the problem it is natural to take 
c = c a at r = a and c = q, at r = b and the zero 
flux of other components dc p /dr = dc p /dr = 0. 
But strictly speaking, the boundary conditions for all 
three constituents in real membranes can not be spec- 
ified basing on a set of biological measurements 
As it can be demonstrated, the specification of one 
or other type of boundaries does not affect cardinally 
the shape of the solution across the whole width of the 
membrane exept thin layers at the edges. In any case 
at chosen way of handling the problem (see the next 
chapter) the set of boundary conditions should be in- 
cluded in the solution by some kind of self-matching 
procedure. 

Integrating the sum of two last equations in (2) yields 



2 Reaction-Diffusion System in Spherical Geom- 
etry 

We start with considering the problem (2)-(3) 
within the spheroidal shell of internal and external 
radii a and b. The consideration of this problem 
within plain and cylindric geometry was performed 
not long ago [ 1 1 and was intended to describe the 
problem of oxygen saturation in the muscular tissue. 
Our choise of the spherical form of the membranes, 
besides purely mathematical interest, is motivated by 
several reasons. Besides the interest in elucidating 
the biological problem of facilitated oxygen transport 
in the (spherically shaped) alveols we refer also to 
the above mentioned diffusion problem in crystalline 
bodies. Usually the mathematical models of impurity 
diffusion use the plain geometry as a tool to repre- 
sent the boundary of a solid body as the platform for 
the diffusion in a bulk. But the surface of a crystal 
is by no means plain, and its inhomogeneities can af- 
fect the diffusion effects (intuitively it is clear from 
the fact that merely the effective diffusive surface is 
bigger than the mathematical surface of the body). 
So, plain model for the boundary is an approxima- 
tion which can be improved. Namely, it is possible 
to model the surface as a highly irregular sequence 
of spheres of different radii, perhaps, forming fractal 
structure. Therefore the consideration of the spheri- 
cal geometry (on a single sphere) is believed to yield 
some improved approximation than just considering 
the straight plane. 

So in the following we speak on the biological prob- 
lem of oxygen transport and consider the spherical 
membrane. We assume the spherical symmetry so 
only one coordinate, namely r resides. As bound- 



c P + c c = const = K 



(4) 



(the conservation of the protein content). If we intro- 
duce the new function 



Y = c c /K 



(5) 



which has a meaning of the ferment saturation func- 
tion the system of equations is cast as: 



D 1 d / 2<9c, _ 
r 2 dr dr 
^ _ Id, 2 dY. 



with 



p = hKc-Y -Kfac + k-) 



(6) 
(V) 

(8) 



Keeping in mind everything said about the bound- 
ary conditions specifications we now try to handle the 
problem setting 



Y 



\a 



Y a Y {b 



C\b = c b 
Y b 



(9) 



where c b and both Y a , Y b should be further deter- 
mined basing upon restrictions imposed by the bio- 
logical sense of the problem. 

Using the boundary conditions add eqs (6-7) and inte- 
grate twice thus obtaining an expression relating c(r) 
and Y(r): 

D{c b -c(r)) + D p K(Y b -Y(r)) = 

-7— • T~~ • [D{c b - c a ) +D p K (Y b - Y a )\ (10) 
br b — a 



3 External Solutions and Facilitated Transport 

Expressing Y from (10) and inserting it into (6) 
and (8) following equation results: 



The variation of the saturation function and the total 
oxygen flux are: 



1 d 2® c \ k 
r 2 dr dr D 
k^D 



D^iz(r 2 ^:) = ~ ■ [Dc b + Y b KD p ] + 



D r , 



hK[Y b - 1 + 



D 



D P K 



Cb\ 



b r ab (k_ [D{Ca _ Cb)+DpK{Ya _ Yb)] _ 
'^[D(c a - c b ) + D p K(Y a - Y b )]\ + ^-c 



br b — a [ D 
k-\ ■ c 



The values Y b , Y a and c b , c a as said above should be 
determined from self-matching conditions. 
The equation (TTTt is extremely difficult (if not impos- 
sible) to solve analytically. Instead of looking for its 
exact solution we notice that <fTTT> can be cast (after 
scaling variables) as the equation in the form: 

eAc = f(c,r), 

and e appears to be a small parameter (~ D/k\). 
This parameter is smaller if the reaction is more in- 
tensive; for biological issues its value varies within 
10~ 4 — 10 -6 (JJ, which allows approximate pertur- 
bative treatment of the problem. Since this param- 
eter enters the equation near the term of the highest 
derivative, the problem appears to be that of the sin- 
gular perturbation theory. Dropping out the e term 
one gets so-called "external solution" to the problem. 
In general this equation of smaller order can't satisfy 
the boundary conditions and need "suturing" with 
the internal (exact) solution around the boundaries. 
The external solutions may suit for practical purposes 
within the bulk of the body, and would yield the de- 
sired solution of (flTt provided that the constants (like 
Y b , Y a , c b etc) are chosen properly, namely by means 
of the self-matching procedure. In this case (fTTl turns 
out to be an ordinary algebraic equation, whose coef- 
ficients depend on the values of Y b , Y a , c b , c a . But in 
fact these latter themselves are solutions to the equa- 
tion at the boundaries. To ensure that this is the case 
one must set following self-matching conditions: 



Y b 



hc b 



k\c h + k. 



Y, 



klC a 



k\c a + k- 



(12) 



(These expressions could be easily obtained from (8) 
setting p = at the boundaries). 



Y n ~ Y h 



kik-(c a - c b ) 



(k\c a + k-){kic b + k-) 



(13) 



and 



(11) 



F{r) 
F d + F f 
D„ ■ K ■ 



-D- 



dc 
) — 

dr 
ab 



k±k- 



dY 

KD >1F S 

[D{c a - c b )+ 

,C a ~ C b ) 



(14) 



{k\c a + k^){k 1 c b + 

which is the sum of two terms, first of which, is 
the ordinary diffusive flux, and the second one, Ff is 
the additional flux caused by the reaction. 
The self-matching conditions yield two relations be- 
tween four values c a , q>, Y a and Y b (or equally s). 
From the point of view of experimentalists these val- 
ues are not treatable in the same fashion. Depending 
on whether we consider the biological problem of fa- 
cilitated oxygen flux or mathematically similar prob- 
lem of impurity diffusion, the set of good-to-operate 
values changes. For the oxygen transfer we can op- 
erate only the value c a within wide range of mag- 
nitudes, and measure s, but this latter only for ex- 
tremely big Ca, where it is known to tend to a constant 
s (~ 0.6 — 0.8) Q. And thus in order to get fluxes 
and concentrations within moderate c a values one re- 
curses to an interpolation which is expected to yield 
semi-quantitative results. Let us express c b from ill 3b : 

i 9 

fc_ic (l - s) 



Cb = T m t- ^ m \ ■ (15) 

k\c a s + fc_(l + s) 

This expression is exact. Hovewer the flux s is itself 
a function of c a to be determined (only its value s at 
c a — > oo is available), so let us go to the semiqualita- 
tive interpolation: 

k-C a (l — s) 

kiC a S + + S) ' 

The expression (TToT i: a) satisfies the limiting condi- 
tion Cb —> at c a = 0; b) at big c a {k\c a » kJ) is 
consistent with dl5t . For intermediate c a this is be- 
lieved to be a plausible interpolation. So from ( TToT) 
the approximate expression for saturation function is 
s = Y a - Y h : 



k\c a s 



k\c a + Ik- 



kic a + k-(l + s)_ 



(17) 



From <fT6b - dT7T >. substituting it into (114b . we get the 

expression for the complete flux through a mem- 
brane. The flux at the external (r = b) layer equals 

F = F d + F f ~ 

-4vr^-^-c a s(A; 1 c a + 2A;_) x (18) 
o o — a 

_/cic a s + k-(l + s) 

{k lCa + k_){k lCa + k_(l + s)) 

and is schematically depicted on Fig. 1 where the or- 
dinary diffusion flux and facilitated one are shown as 
functions of c a . 



F 




Figure 1 : Diffusion and reaction facilitated oxygen 
flux as function of c a (arbitrary units). 



It is interesting to note different character of two 
fluxes from the latter formula and from Fig. 1 : at big 
c a the diffusion flux grows linearly, but the flux com- 
ponent due to selfconjugation with reaction tends to 
a constant value. For both biological problem, and 
for the problem of vacansion diffusion enhancement 
this is clear intuitively: since the amount of carrier is 
limited at big values of c a the substrate simply gets 
saturated. As to the moderate or small c a the relation 
of two fluxes is 

F f _ D p Kki 

T d Dk^ ■ 

The value D S> D p (since D stands for "fast 
species"), but under condition k\ k- which is true 
far from chemical equilibrium, and if it is the case, 
the facilitated transport dominates over diffusion. 




Figure 2: Space concentration of the free component 
c(r) and total diffusion component c(r) + c c (r) . 

In the solid body physics the role of the reversible re- 
action is undertaken by the process of reversible tran- 
sitions of impurity atoms between interstitial and sub- 
stitution positions. The process of this interchange 
can be represented as the chemical reaction with va- 
cancy: 

/ + y — S (19) 

(I and S standing for interstitial and substitutional 
impurities, V for vacancies). Thus vacansies can be 
formally understood as "impurity carriers" like Hb 
complexes. Another mechanism is so called Watkins 
mechanism which also involves the interchange of 
fast and slow species, but by the intermittency of 
eigen interstitial atoms of the matrix. It is known 
(e.g., [2 1) that such reversible reactions are crucial 
for understanding the impurity redistribution in the 
depth of a crystal. In order to handle analytically 
the impurity profile, the authors of [2| used similar 
assumptions of rapid chemical equilibrium establish- 
ing compared to the ordinary diffusion. In conclusion 
on Fig.2 we show the approximate solution for space 
oxygen concentration in a spherical slab (between in- 
ner a and outer b boundary). Shown are both diffu- 
sion c(r) component and the value c + c c , that is total 
oxygen concentration. 
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